args <- commandArgs(TRUE)
setwd(args[1])
var<-read.table("inputfile.txt",sep="\t")
library(affy)
vec1<-as.matrix(aggregate(var[,2], list(var[,1]), tukey.biweight))
#colnames(vec1)<-c("Gene","expr")
vec1<-cbind(vec1,rep(0.05,nrow(vec1)),rep("T",nrow(vec1)))
colnames(vec1)<-c("SEGID","EXPR_MEAS","CONF","CALL")
#m<-noquote(c("C","A","K"))
#m1<-paste(m,collapse="\t")
#write.table(m1,"sample.t",col.names=FALSE,row.names=FALSE,sep="\t",quote=FALSE)
write.table(vec1,"outfile.txt",sep="\t",row.names=FALSE,col.names=TRUE,quote=FALSE)
q()
n
